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Abstract. The relevant energy ranges for stellar nuclear reactions are introduced. Low-energy com- 
pound and direct reactions are discussed. Stellar modifications of the cross sections are presented. 
Implications for experiments are outlined. 
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INTRODUCTION 
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Astrophysical reaction rates describe the number of reactions in a given volume of the 
stellar plasma and are the essential ingredient for calculations following nucleosynthe- 
sis and energy generation in astrophysical environments. They are derived by folding 
experimental or theoretical reaction cross sections with the energy distribution appro- 
priate for the interacting plasma constituents. For nucleons and nuclei in stars this is the 
Maxwell-Boltzmann distribution (MBD). A direct determination of the cross sections 
faces several problems. The majority of nuclei participating in explosive nucleosynthe- 
sis is short-lived and thus the reactions cannot (yet) be studied in the laboratory. Ther- 
monuclear reactions on stable and unstable targets proceed at low, often subCoulomb, 
energies, resulting in tiny reaction cross sections which are further obstacles to experi- 
ments. On the theoretical side, several reaction mechanisms may contribute, each with 
its own peculiar uncertainties. This makes the determination of reaction rates at astro- 
physically relevant energies a complex puzzle which currently can only be tackled by 
using several different approaches. In the following I briefly outline a few aspects im- 
portant for all investigations. 



COMPOUND REACTIONS 

The majority of astrophysically relevant reactions proceeds through a compound nu- 
cleus (CN). We can distinguish between compound reactions exhibiting pronounced 
resonance features in their excitation functions and such without those features. The nu- 
clear level density (NLD) at the compound formation energy (i.e. projectile energy plus 
separation energy of the projectile in the CN) is the distinguishing factor. A high NLD 
results in a large number of tightly overlapping resonances which cannot be resolved 
whereas a low NLD leads to more or less isolated resonance structures. Systems with a 
high NLD can be treated in a statistical approach, assuming an average over the NLD. 



This is called the Hauser-Feshbach model (HFM). Isolated resonances (and their inter- 
ference) have to be treated in phenomenological approaches such as the Breit-Wigner 
formula (BWF) or R-matrix fits. 



Transmission coefficients and widths 

The reaction cross section of the CN reaction a+A— > b + B is defined by summing 
over strength functions, widths or transmission coefficients 

n Mot 

where the sum is over m individual resonances (BWF, n = 1 . . . m) or over all spin and 
parity pairs in the compound system at its formation energy (HFM). In the BWF the 
Ts in the numerator are the partial widths of the resonances in the entrance and exit 
channels. In the HFM these are averaged widths T = (T) which are related to strength 
functions S? — p(T) and transmission coefficients T = 2np(T) = 2% 5^ '. The density of 
levels with given spin J and parity % in the compound system at the formation energy is 
denoted by p . In both cases, r tot is the total width computed from the sum of all partial or 
averaged widths for all open reaction channels (including re-emission of the projectile). 
It is noteworthy that Tf, (as well as r tot ) involve sums over transitions to all possible 
final states in the given channel whereas Y' a only includes the transitions from the target 
ground state (g.s.) when laboratory cross sections are to be calculated. Additionally, all 
widths for a given transition include sums over the allowed partial waves. 

Equation (1) and the (averaged) widths entering are essential to understanding a 
number of intricacies connected to the theoretical and experimental determination of 
cross sections for astrophysical reaction rates. This is outlined in a few examples in the 
following chapters. 



Relevant energy ranges 

Nuclear reactions in stellar plasmas proceed at interaction energies which are low 
by nuclear physics standards, even in explosive nucleosynthesis. The energy range 
contributing most to a stellar rate r at given stellar temperature T is found by inspection 
of the integral arising from folding the reaction cross sections with the MBD 
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Although the integration formally runs from Zero to Infinity most of the contributions 
to the integral come from a narrowly defined energy region Eq ± (AE)/2. For charged 
particles this relevant energy region is usually called the Gamow window. The relevant 
energy region is determined by the energy dependences of the cross sections and the 
MBD. Frequently used are simple approximation formulas for the Gamow window (e.g., 



as given in [1, 2]). These are based on the implicit assumption that Y' a <C Tj,. It has been 
pointed out that this assumption is not always fulfilled, especially for intermediate and 
heavy target nuclei [2, 3,4]. Using Eq. (1) it is easy to see why. If the largest width in 
the numerator is also dominating r to t it will cancel with the denominator and the energy 
dependence will be solely given by the smaller width in the numerator. If this happens 
to be the entrance (projectile) width, the application of the standard approximation 
formulas is justified. However, in many cases it is not. 

A straightforward example is the one of capture. At sizeable projectile energies, the 
7 width will always be smaller than the projectile width. This will certainly apply to 
neutron capture. Fortunately, the energy dependence of s-, p-, and d-waves is much 
weaker than the exponential behavior of the MBD and in consequence the relevant 
energy window is given by the peak of the MBD. The fact that the y width is the smallest 
also holds for charged-particle capture, unless with light targets or at low interaction 
energies. In the latter case, the projectile width can be suppressed by the Coulomb 
barrier and may become smaller than the photon width. This depends on the structure 
of the nucleus and the reaction energetics and has to be scrutinized separately for each 
reaction. Another complication is to have different Coulomb barriers in the entrance 
and the exit channel. Which width is the smallest will then sensitively depend on the 
projectile energy and the reaction Q value. 

A detailed investigation of these issues and newly derived ranges of astrophysically 
relevant energies can be found in [4] . Shifts by several MeV of the energy windows have 
been found. It is strongly recommended to consider these energy windows instead of the 
ones calculated with the standard approximation formulas, especially when devising 
experiments of astrophysical interest. 

EXPERIMENTAL CONSIDERATIONS 

Determination of astrophysically relevant cross sections 

Equation (1) is also central to the question of how much astrophysically relevant 
information can be extracted from a measurement of reaction cross sections. Foremostly 
it has to be attempted to fully cover the relevant energy window with measurements. 
The energy windows given in [4] are defined in such a way that knowledge of the cross 
sections across the full energy range will determine the reaction rate integral (Eq. 2) to 
10%. For higher accuracy, measurements have to extend to higher and lower energies. 

Frequently, it is not possible to fully cover the relevant energy range. Especially for 
charged-particle reactions at intermediate and heavy targets, often it is not even possible 
to reach the Gamow window because subCoulomb reactions exhibit tiny cross sections. 
In this case measurements at higher energy are extrapolated towards lower energy or 
are compared to calculations to check the theoretical models. Caution is advised in 
these attempts because the energy dependence of the cross section and the sensitivity 
of the calculation to certain nuclear inputs can change depending on the energy under 
investigation. Again, Eq. (1) provides guidance. Energy dependence and sensitivity are 
determined by the smallest width in the numerator. (Additional complications arise when 
both widths are comparable or another channel is considerably contributing to r to t.) 
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FIGURE 1. Sensitivities s for 96 Ru(p,y) 97 Rh 



Since the sizes of the widths relative to each other change with energy, the smallest 
width may be found in another channel at low energy than at high energy. 

A helpful visual aid to estimate the relative importance of the different channels is a 
sensitivity plot. We define the sensitivity s as a measure of a change in the cross section 
fa = Chew/Cold as the result of a change in a width by the factor f m , with 5 = when 
no change occurs and s = 1 when the cross section changes by the same factor as the 
width: 
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Plotting s as a function of the cm. energy yields a plot like the example shown in Fig. 
1. Using this example for investigating the reaction 96 Ru(p,y) 97 Rh, we find from [4] 
that the Gamow window is 1.63 < E < 3.42 MeV for the typical p-process temperature 
T = 2.5 GK. It can clearly be seen in Fig. 1 that the sensitivities are very different 
at lower and higher energies. For example, a measurement closely below the neutron 
threshold would be in a region where s is largest for the y width but smallest for the 
proton width, just the opposite of what is found in the astrophysically relevant energy 
region. Above the neutron threshold the situation is even more complicated because there 
is additional sensitivity to the neutron width (dominating r tot ), although not as large as 
to Ty. If any discrepancy between measured and predicted cross sections was found, it 
would be hard to disentangle the different contributions. In any case, no information on 
the astrophysically important proton width could be extracted from a measurement at 
the higher energies shown in Fig. 1. 

From the general considerations above it also follows that it is advantageous to use 
reactions with neutrons in one channel for investigating the sensitivity to the charged- 



TABLE 1. Stellar enhancement factors (SEF) for selected (y,n) reactions from NON- 
SMOKER calculations [7] 

Target: SEF 

186 W:400 185 Re: 1300 187 Re: 1200 190 Pt: 5500 192 Pt: 3300 198 Pt: 310 
197 Au: 1100 196 Hg: 1700 198 Hg: 750 204 Hg: 43 204 Pb: 160 



particle optical potential, i.e., using (a,n), (n,a), (p,n), or (n,p) reactions. Except within 
a few keV above the neutron threshold, the neutron width will always be much larger 
than the charged-particle width, even at higher than astrophysical energies. Therefore it 
will cancel with the denominator in Eq. (1) and leave the pure energy dependence of the 
charged-particle width. On the other hand, information on the neutron potential is hard 
to get from reactions. But this also implies that the sensitivity of astrophysical rates to 
the neutron optical potential is not high. 

On a side note, if averaged widths (I") or strength functions S? can be determined 
experimentally, it is more useful to obtain the latter. They are directly proportional to the 
transmission coefficients T calculated in the theoretical reaction models by solving the 
Schrodinger equation, without an additional dependence on the compound NLD. 



Stellar Effects 

Astrophysical investigations require the use of stellar rates accounting for thermal 
population of excited states in the target. Only stellar rates allow to convert the rate 
for the forward reaction to the one for the reverse reaction by a simple relation de- 
rived from detailed balance. Stellar rates r* are calculated using stellar cross sections 
<?* = (L M (2 M + l)exp(- J E M /(fcr)(7 M )/(L M (27 /J + l)exp(-£ M /(fcr)) (where \i indi- 
cates the excited target states) instead of laboratory cross sections o — On=o- It can 
be shown [5, 6] that this leads to an equation identical to the one in Eq. (2) but with a 
replaced by an effective cross section a eff , including sums over transitions to all ener- 
getically accessible states in the exit and entrance channel. This means that Y' a in Eq. 
(1), including only transitions to the target g.s., is replaced by a T a similar to Tf,. 

The stellar enhancement factor (SEF) S = o*/o is useful to assess the impact of the 
excited states. It is also useful to see how far the measured cross section is from the 
desired stellar cross section. A measurement should always be made in the direction of 
the SEF as close as possible to Unity, the stellar cross section (or rate) of the inverse 
reaction can then be calculated from detailed balance. From the definition of c eff and 
Eq. (1) it is easy to see that this will usually be the direction of positive reaction Q value 
as in this case fewer transitions will contribute in the entrance channel than in the exit 
channel. Photodisintegration reactions are the most extreme cases because they exhibit 
the largest difference in the number of contributing transitions between the entrance 
and exit channel. Table 1 gives an example of the SEF appearing in such reactions 
[7, 8]. Clearly, transitions from and to the g.s. (which are measured in the laboratory) 
contribute only a small fraction of the stellar cross section and it is better to measure 
capture reactions in order to be as close as possible to the stellar value. This is because 



the NLD is strongly increasing with excitation energy whereas the 7 strength function 
is decreasing with decreasing 7 energy. This leads to a maximal contribution to the 7 
strength of 7 transitions 3—4 MeV below the compound formation energy [9, 10]. 

There are some exceptions to the rule of measuring in the direction of positive Q 
value, though [11, 12]. The range of transition energies in each channel varies from Zero 
to the particle separation energy. Transitions with a small relative energy are strongly 
suppressed by the presence of a barrier, e.g., the Coulomb barrier or the centrifugal 
barrier. If entrance and exit channel exhibit very different barriers, an endothermic 
reaction may show a smaller SEF than its inverse reaction although in principle the 
number of possible transitions is larger. Most of these transitions, however, may be 
suppressed by the barrier. The occurrence of this phenomenon depends on the barriers in 
the entrance and exit channels, and on the reaction Q value. A very detailed explanation 
of the effect, along with a list of affected reactions, can be found in [12]. 

REACTION MECHANISMS 

Modification of the statistical model 

The prediction of astrophysical reaction rates is difficult not just because the required 
nuclear properties are not fully known but also because several reaction mechanisms 
may contribute. Most astrophysically relevant reactions can be described in the HFM. 
Regular HFM calculations assume a compound formation probability independent of 
the compound NLD at the compound formation energy. Therefore the sum in Eq. (1) 
runs over all J n pairs (a high-spin cutoff is introduced in practical application of the 
model because spin values far removed from the spins appearing in the initial and final 
nuclei do not contribute to the widths). The availability of compound states and doorway 
states defines the applicability of the HFM [13]. Relying on an average over resonances, 
the HFM is not applicable with a low NLD at compound formation. Single resonances 
and direct reactions will contribute then. On average the HFM will then overpredict the 
resonant cross section (unless single resonances dominate) because it will overestimate 
the compound formation probability. This can be treated by introducing a modification 
of the formation cross section which includes the compound NLD dependence. The 
summands of Eq. (1) will then be weighted according to the available number of states 
with the given J n . (Formally this is the same as assuming J n dependent potentials for 
particle channels.) For an implementation of a parity dependence this was discussed 
in [14]. This was also used in the HFM code NON-SMOKER WEB since version 4.0w 
[15]. Additionally, the option of weighting the HFM cross section by a NLD dependent 
function was offered. The new code SMARAGD [16, 17, 18] has an improved version 
of this as default and uses J K dependent weighting of the summands, thus implicitly 
accounting for a low NLD at the compound formation energy. Preliminary results with 
this modification are shown in Fig. 2. 
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FIGURE 2. Neutron capture cross sections of even Sn isotopes in an averaged DC model with energy- 
dependent spectroscopic factor (DC), a modified HFM (mod HF), and a standard HFM (HF) (preliminary 
results) 



Low-energy direct reactions 

Due to the low particle separation energies and/or NLDs encountered in very neutron- 
(or proton-) rich nuclei, direct reactions will also become important, even at very low 
energies (in the keV region) [13, 19, 20]. (The direct-semidirect mechanism at higher 
energies has been found to be of minor relevance for astrophysics [21] .) Mainly relevant 
for astrophysical applications is direct capture (DC). The relative importance of DC in- 
creases with decreasing neutron separation energy S n for (n,y) on neutron-rich nuclides 
[19] and with decreasing proton separation energy S p for (p,y) on proton-rich targets. 
This is because the number of final levels for y transitions from resonant capture be- 
comes lower [9]. The Sn isotopes have implicitly a very low NLD and consequently the 
highest relative importance of DC. For this reason, they have become the focus of many 
experimental and theoretical investigations [9, 10, 14, 19]. 

The two major problems in the prediction of DC lie in the determination of the 
final states (i.e., the low-lying compound states) and their spectroscopic factors far off 
stability. Nuclear structure models cannot predict these states with sufficient accuracy 
and predicted DC cross sections differ by several orders of magnitude [22] . Furthermore, 
most DC predictions far off stability have utilized only constant spectroscopic factors so 
far [23], only few have employed energy-dependent ones [24, 25, 26]. To circumvent 
the problem of the exact prediction of final states it has been suggested [24, 25, 23] 
to employ averaged properties, i.e. to replace the sum over discrete final states by an 
integration over a level density, similar to the HFM. (The determination of the NLD 



poses the same problem as in the case of compound reactions.) 

The code SMARAGD also includes a global DC treatment using an averaged DC 
model and energy-dependent spectroscopic factors based on BCS and Lipkin-Nogami 
occupation probabilities [21, 19, 17, 19]. This average DC approach aims at providing 
robust predictions despite of considerable differences between microscopic predictions 
[22]. Preliminary results for this DC treatment are shown in Fig. 2. 

The final rate (or cross section) is the sum of the modified HFM value and the DC 
one. Interestingly, for the Sn isotopes shown here (except for TV = 92) this sum is 
approximated by the unmodified HFM result within a factor of 10. This shows that it 
seems justified to use unmodified HFM rates as crude estimate of the total rates for 
exotic nuclei. 

This work was supported by the Swiss NSF, grant 200020-122287. 
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